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Abstract 

The influence of compressibility and helicity on the stability of the scaling regimes of a passive 
scalar advected by a Gaussian velocity fleld with finite correlation time is investigated by the 
field theoretic renormalization group within two-loop approximation. The influence of helicity and 
compressibility on the scaling regimes is discussed as a function of the exponents e and 77, where 
e characterizes the energy spectrum of the velocity field in the inertial range E oc fc^~^^, and rj is 
related to the correlation time at the wave number k which is scaled as /c~^+''. The restrictions 
given by nonzero compressibility and helicity on the regions with stable infrared fixed points which 
correspond to the stable infrared scaling regimes are discussed. A special attention is paid to the 
case of so-called frozen velocity field when the velocity correlator is time independent. In this 
case, explicit inequalities which must be fulfilled in the plane e — rj are determined within two-loop 
approximation. 

Introduction 

One of the main problems in the modern theory of fully developed turbulence is to verify the validity 
of the basic principles of Kolmogorov-Obukhov (KO) phenomenological theory and their consequences 
within the framework of a microscopic model On the other hand, recent experimental, numerical 

and theoretical studies signify the existence of deviations from the well-known Kolmogorov scaling 
behavior. The scaling behavior of the velocity fluctuations with exponents, which values are different 
from Kolmogorov ones, is known as anomalous and is associated with intermittency phenomenon [2]. 
Even thought the understanding of the intermittency and anomalous scaling within the theoretical 
description of the fluid turbulence on basis of the "first principles", i.e., on the stochastic Navier-Stokes 
equation, still remains an open problem, considerable progress has been achieved in the studies of the 
simplified model systems which share some important properties of the real turbulence. 

The crucial role in these studies is played by models of advected passive scalar field . Maybe the 
most known model of this type is a simple model of a passive scalar quantity advected by a random 
Gaussian velocity field, white in time and self-similar in space, the so-called Kraichnan's rapid-change 
model It was shown by both natural and numerical experimental investigations that the deviations 
from the predictions of the classical KO phenomenological theory is even more strongly displayed for 
a passively advected scalar field than for the velocity field itself (see, e.g., jS] and references cited 
therein). At the same time, the problem of passive advection is much more easier to be consider from 
theoretical point of view. There, for the first time, the anomalous scaling was established on the basis 
of a microscopic model ||6j, and corresponding anomalous exponents was calculated within controlled 
approximations (see review [SI and references therein). 
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In paper [Jj the field theoretic renormalization group (RG) and operator-product expansion (OPE) 
were used in the systematic investigation of the rapid-change model. It was shown that within the 
field theoretic approach the anomalous scaling is related to the very existence of so-called " dangerous" 
composite operators with negative critical dimensions in OPE (see, e.g., [SUHl for details). 

Afterwards, various generalized descendants of the Kraichnan model, namely, models with inclusion 
of large and small scale anisotropy ^U], compressibility and finite correlation time of the velocity 
field jl21 I13j were studied by the field theoretic approach. General conclusion is: the anomalous 
scaling, which is the most important feature of the Kraichnan rapid change model, remains valid for 
all generalized models. 

In paper the problem of a passive scalar advected by the Gaussian self-similar velocity field 
with finite correlation time [T^ was studied by the field theoretic RG method. There, the systematic 
study of the possible scaling regimes and anomalous behavior was present at one-loop level. The 
two-loop corrections to the anomalous exponents were obtained in \'> . In paper jl3j the influence of 
compressibility on the problem studied in ^21 was analyzed. The effects of the presence of helicity 
in the system on the scaling regimes and anomalous dimensions were studied in |161 117j within of 
two-loop approximation. It was shown that although the separate composite operators which define 
anomalous dimensions strongly depend on the helicity parameter the resulting two-loop contributions 
to the critical dimensions of the structure functions are independent of helicity. This rather intriguing 
result can have at least two independent explanations. First, it can be simply a two-loop result which 
will be changed when one will study three-loop approximation. Second, more interesting conclusion 
(but it is only a speculation for now) is that this situation will be held at each level of the perturbation 
theory while the assumption of incompressibility or isotropy will be supposed. From this point of view, 
the investigation of the compressible system together with assumption of helicity of the system within 
two-loop approximation can give an interesting answer. For this purpose in jl8j the influence of 
compressibility of the model on the stability of the scaling regimes was studied and the restrictions 
on the parameter space were analyzed. In what follows, we shall continue these studies and our aim 
will be to find combine effects of helicity and compressibility on the stability of the scaling regimes in 
two-loop approximation. It can be consider as the starting point for the subsequent investigation of 
the anomalous scaling of the correlator or structure functions of a passive scalar. 



Formulation of the model 

The advection of a passive scalar field is described by equation 

dt9 + (v • d)e = u^ve + f , (1) 

where 0{x) = 0{t,x.) is a passive scalar field, z^o is the molecular diffusivity coefficient, v(x) is com- 
pressible velocity field, and = f^{x) is a Gaussian random noise with zero mean and correlation 
function 

{f{x)f{x')) = 5{t - t')C{r/L), r = X - x', (2) 

where parentheses (...) denote average over corresponding statistical ensemble. The noise maintains 
the steady-state of the system but its concrete form is not essential. L denotes an integral scale related 
to the stirring. 

In real problems the velocity field v(x) satisfies stochastic Navier-Stokes equation but we shall 
suppose that the velocity field obeys a Gaussian distribution with zero mean and correlator 

{v,{x)v,{x')) =J^J ^i?.,(k)D„(^, A:)e-(*-*')+^'^-(--') . (3) 

Here -Rjj(k) = Pij{\i.) + (5ij(k) -|- i7jj(k), where -Pij(k) = 5ij — kikj/k^ is the non-helical transverse 
projector, Qij(k) = akikj/k'^ is longitudinal projector with the compressibility parameter a > 0, 
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and Hij{h) = ipeijiki/k is the helical transverse projector with helicity parameter p E (0,1) {siji is 
Levi-Civita's completely antisymmetric tensor of rank 3, and k = |k|). d is the dimensionality of the 
X space which must be taken d = 3 in the helical case {p > 0). The function Dy is taken in the form 

Dy{u,k)= (4) 

where 50 > and the exponent e describe the equal-time velocity correlator (the energy spectrum). 
On the other hand, the constant uq and the exponent r] are related to the frequency lo, characteristic 
of the mode k. In our model the exponents e and r] play the role of the RG expansion parameters. The 
model contains two special cases. In the limit — > 00 and Qq/uq = const, one obtains the so-called 
rapid-change model, and the limit ^ and go/uQ = const, corresponds to the case of a "frozen" 
velocity field [l2|IlB]. 

The stochastic problem (EO)-® is equivalent to the field theoretic model of the set of fields $ = 
{0,0', v} (see, e.g., (Hj) with action functional 

S{^)= - a!tid'^xidt2<i'^X2 ^;i(ti,xi)[L»j;.(ti,xi;f2,X2)]-^i;j(t2,X2) 

+ J dt d'^x 9' [-dtO - VidiO + uqAO] , (5) 

where the unimportant term related to the noise (j2I) is omitted, 9' is an auxiliary scalar field, and 
summations are implied over the vector indices. The second line in represents the Martin-Siggia- 
Rose action for the stochastic problem at fixed velocity field v, and the first line describes the 
Gaussian averaging over v defined by the correlator in (jS)) and @. 

The functional formulation © means that statistical averages of random quantities in the stochas- 
tic problem defined by and © corresponds to functional averages with the weight S{^). 

UV renormalization, fixed points and scaling regimes 

Using the standard analysis of quantum field theory one can find that the UV divergences are only 
present in the one-particle-irreducible Green functions {9'9)i-ir and (0'0v)i_jj.. The renormalized 
action functional has the following form 

S{^)= - ^y"dtid%dt2d'^X2t;i(ti,xi)[i?,];.(ti,xi;t2,X2)]-^^;j(t2,X2) (6) 
+ J dtd'^^ 9' [-dt9 - Z2Vidi9 + vZiM] , 

The dimensionless parameters 5, u and u are the renormalized counterparts of the bare parameters 
(denoted by 0). Zi and Z2 are the renormalization constants of our model. We shall calculate them in 
2-loop approximation in the minimal subtraction scheme. Our model is multiplicatively renormalizable 
which is represented by the following relations between renormalized and bare parameters (see, e.g., 

i/Q = vZy, go = gp^'^^'^Zg, no = up'^Z^, v Z^v. (7) 

where 

Zu = Zi, Zu = Zi ^, Zg = Z\Zy Zy = Z2. (8) 

and p is the renormalization mass. One and two-loop Feynman diagrams which contribute to the 
renormalization constants Zi and Z2 can be found in They have the following form 



Zi = ^Ai + ^(-Bi + cX Z2 = ^A2 + ^(-B2 + C2], (9) 
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where the one-loop contributions are (in the MS-scheme) 

^ Sd 1 {l + u){d-l + a) -2a ^ Sd a , . 

^~ (27r)'^ 4u(l + n)2 d ' ^ ~ (27r)'^ 4«(1 + «)2 ' ^ ' 

where Sd is d dimensional sphere given by the expression Sd = 27r'^/2/r(d/2). Two-loop contributions 
will not be shown explicitly here (the explicit form of the corresponding expressions in the non-helical 
case can be found in ^H], and the explicit form of the helical two-loop contribution to the Zi can be 
found in JH])- 

From the renormalization constants @ one can define the corresponding anomalous dimensions 
7i = fidfj^ In Zi and /3 functions for all invariant charges {X = g,u): (3x = fJ-d^X. In our case we have 

71 = /uS^lnZi = -2igAi + 2g^Ci), 72 = fid^lnZ2 = -2{gA2 + 2g^C2), (11) 

and 

l3g = fJ'd^tg = g{-2e - T] + 3-fi - 2^2), Pu = fJ'd^u = u{-ri + ji). (12) 

Now one has all necessary tools to begin with analysis of possible scaling regimes of the model. It 
is well known that possible scaling regimes of a renormalizable model are associated with the infrared 
(IR) stable fixed points of the corresponding RG equations. The fixed points are determined from the 
requirement that all /3 functions are vanish: /3x(-'^*) = 0, where X = g,u and X^ are coordinates of 
the corresponding fixed point. The type of the fixed point is given by the matrix of the first derivatives 
Qij = d(3i/dXj. The fixed point is IR stable if all the eigenvalues of the matrix are positive (precisely 
their real parts). In our model we have five types of possible scaling regimes related to fixed points of 
the model. 

First of all, we shall investigate the rapid-change limit: u 00. In this regime, it is necessary 
to make transformation to new variables, namely, w = 1/u, and g' = g/u^, with the corresponding 
changes in the (3 functions: 

/Jg. =g'(?7-2e + 7i- 272), (3^ = w{r] - li)- (13) 

In this case we have two fixed points which we denote as I and II. First of them is trivial one, namely 

Fixed Point I: = l/u^, = g^ = 
with 71 = 0. The matrix 0, is diagonal with the elements (eigenvalues): 0,i = rj, 0,2 = V ~ 2e. Thus, 
the corresponding scaling regime is IR stable if 77 > and, at the same time, 77 > 2e. 

The second one is defined as follows 

Fixed Point II: tu* = l/w* = 0, and 

, ]\ (2g-r?), (14) 
d — 1 + a 

where we denote g' = g'Sd/ (27:)'^. The anomalous dimension 71 at the fixed point is 7J = 2e — rj. These 
are exact one-loop results. It is the consequence of the non-existence of the higher-loop corrections. 
Corresponding matrix of the first derivatives is triangular with diagonal elements (eigenvalues) f^i = 
2{rj — e) and ^2 = 2e — rj. Thus, the conditions g[ > and ^}i^2 > for the IR stable regime lead to 
the inequalities: rj > e and r/ < 2e. 

The second nontrivial limit of our general model is so-called frozen limit given by n ^ 0. In 
this case, it is again necessary to make transformation to new variables, namely, g" = g/u and u is 
unchanged. It leads to the changes in the [3 functions which are now define as follows {I3u is the same 
as in the general case) 

=/(-2e + 27i- 272), Pu = u[-rj + ^i). (15) 

In this case we have again two fixed points which we denote as III and IV. First of them is trivial one, 
namely 
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Fixed Point III: = g'l = 0. 

The corresponding matrix Q is diagonal with the elements: = — 2e and ^2 = — ??• It means that 
the corresponding scahng regime is IR stable if e < and t] < 0. 

The second one is defined as follows 

Fixed Point IV: = 0, and 



n" = ! ^20 2 ^g^ 

y* <)( All _ All \ n( All _ All ^3 ' ^ ' 



where A'-^ and C'/q for i = 1,2 are expressions from Q taken in the given hmit and again g" = 
g" Sd/{27r)'^ . In this case we have the relation = 72 between anomalous dimensions. The matrix 
of the first derivatives is triangular with eigenvalues (diagonal elements): 

In the helical case, when helicity parameter p > 0, one must work with d = 3, and we obtain 

-9': = 3^+^'(^-r-|-v), (18) 

= e-r/-^ + ^(-64 + 56a + 97rV), = e (^2 + e (^-1 + '^-a + ^tt' p^^^^ (19) 
therefore, the conditions to have g'l > and Qi^2 > are 

e>0, e(^-l + Ia + ^7rV) > -2, r? < e - ^ + ^(-64 + 56q + 97rV)- (20) 

These inequalities define the region in the parameter space for which the scaling regime is IR stable. 

The last but the most interesting scaling regime is obtained when one assume that < u^, < 00. 
Let us briefly discuss this case. The corresponding IR fixed point will be denoted as V. The coordinates 
of the fixed point is now defined by the requirement of vanishing of the /3 functions which are given 
in ()12|). The fixed point value for g is given as 

Fixed Point V: > 

^^^^ ^- = -2(1^)- i^x;^^'- P') 

where the functions Ai, A2,Ci, and C2 are given in ®, and where the parameter u is taken at its 
fixed point value which is defined implicitly by the equation 

-V + 7*iM = 0- (22) 

The analysis of the problem is relatively simple in one-loop level, where explicit expressions can be 
found p^. On the other hand, as was briefly discussed in JH], the situation is essentially more compli- 
cated when we are working in two-loop approximation. Detail analysis shows that the investigation of 
the IR stability of the fixed point in the general case of the present model has to be done individually 
for concrete situation which is rather cumbersome and it must be done in a separate work. 



Conclusions 

In present paper we have studied the influence of compressibility and helicity of the system on the 
possible IR scaling regimes of the model of a passive scalar advected by a Gaussian velocity field with 
finite time correlations by means of the field theoretic RG technique. The dependence of the fixed 
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points, which are directly related to the existence of possible IR scaling regimes, on the parameters of 
compressibility and helicity as well as their IR stability is discussed. The explicit inequalities, which 
define the stable IR scaling regimes, are found in the case of the frozen limit of the model. The most 
general case with finite time correlations of the velocity field is more complicated within two-loop 
approximation and has to be consider in detail once more. 
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